Rapid Domain Wall Motion in Permalloy Nanowires Excited by Spin-Polarized 
Current Applied Perpendicular to the Nanowire 
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We study domain wall (DW) dynamics in permalloy nanowires excited by alternating spin- 
polarized current applied perpendicular to the nanowire. Spin torque ferromagnetic resonance 
measurements reveal that DW oscillations at a pinning site in the nanowire can be excited with 
velocities as high as 800 m/s at current densities below 10^ A/cm'^. 

PACS numbers: 75.60.-d, 75.70.-i, 72.25.-b 



Domains in ferromagnetic nanostructures are candi- 
dates for information storage [H Q and processing 
To be technologically useful, magnetic domains must be 
easily manipulated. One way this can be accomplished 
is by moving the magnetic domain wall (DW) separating 
two oppositely magnetized domains. Since manipulation 
of individual DWs in nanostructures with magnetic field 
[1, 0] is technically challenging, DW motion induced by 
spin torque (ST) from spin-polarized current has emerged 
as a promising alternative 16| . The effects of ST on DW 
arc most readily studied in the ferromagnetic nanowire 
geometry with electric current applied along the nanowire 
[THlSf. However, spin torque exerted on a DW in this ge- 
ometry is small because the angle between the current 
polarization direction and magnetization is small every- 
where in the DW. As a result, high current density (>10^ 
A/ cup) is required to achieve DW velocities in the tech- 
nologically useful range of '-^lO^ m/s 

High DW velocities (180 m/s) at relatively small cur- 
rent densities (< 5x10^ A/cm^) have been observed in 
spin valves with current flowing in the plane of the mag- 
netic layers This high efficiency of the current- 
induced DW motion Q was attributed to transfer of an- 
gular momentum from the fixed layer to the DW in the 
free layer of the spin valve mediated by the component of 
spin current perpendicular to the layers. Khvalkovskiy et 
al. [3| also numerically studied DW motion in spin valve 
nanowires with current flowing perpendicular to the spin 
valve layers (CPP geometry). This work predicted that 
DW velocities of ^10^ m/s can be achieved at current 
densities of ^10^ A/cm^ in the CPP geometry [l^. 

In this Letter we make measurements of CPP DW 
dynamics in the permalloy (Py=Ni84Fci6) layer of 
CosoFeso/ Cu/ Py spin valves patterned into 5 /im long 
and 90 nm wide nanowires. The nanowires are defined 
on metallic Cu/Ta/Ru films serving as bottom leads, and 
500 nm wide Ta/Au top leads are used to apply current 
perpendicular to the spin valve layers as shown in Fig. 1. 
The device in Fig. 1 is made in a multi-step nanofabri- 
cation process starting from a Ta(5)/ Cu(30)/ Ta(3)/ 
Cu(30)/ Ta(5)/ Ru(lG)/ Cu(3)/ Co5oFe5o(7)/ Cu(5)/ 
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FIG. 1: (a) Schematic of a CosoFeso/Cu/Py spin valve in the 
shape of a 90 nm wide nanowire with a 500 nm wide top lead. 
A domain wall in the Py layer is trapped at a random pinning 
site under the top lead while the CosoFeso layer is uniformly 
magnetized, (b) Scanning electron microscopy image of the 
device, showing the spin valve nanowire, the top lead and the 
bottom lead, (c) Transverse domain wall profile in the Py 
nanowire obtained from micromagnetic simulations. 



Py(3)/ Cu(5)/ Ru(2.5)/ Ta(2.5) multilayer (thicknesses 
in nanometers). The multilayer is deposited onto a ther- 
mally oxidized silicon substrate by magnetron sputtering 
and annealed at 225 °C for 2 hours. The Py layer magne- 
tization measured with vibrating sample magnetometry 
{Ms= 430 emu/cm-^) is reduced compared to the bulk 
value due to magnetically dead layers at the Cu/Py in- 
terface and interdiffusion of Cu and Py [iot . 

Fig. 2(a) shows the resistance of the device as a func- 
tion of magnetic field applied parallel to the nanowire. 
The giant magnetoresistance (GMR) data in Fig. 2(a,b) 
reveal that the Py layer undergoes reversal in the -1-200 to 
-300 Oe field range whereas the CosoFeso layer switches 
at -900 Oe. Fig. 2(b) shows a major GMR hysteresis loop 
of the Py layer (±600 Oe). An intermediate resistance 
(IR) state in this loop arises from DW trapping near the 
middle of the top lead. The DW is trapped in the same 
state in a minor hysteresis loop (±200 Oe) in Fig. 2(c). 
All measurements reported in this Letter are made at T 



— 80K in order to prevent thermally-activated dcpinning 
of the DW from the IR state observed at 300K. We ex- 
pect the DW velocities at T = 300K to be similar to those 
measured at T = 80K [2lj. For the thin- film nanowire 
geometry of the the Py layer, only a transverse DW is 
expected to be stable [9[. We confirm this th roug h mi- 
cromagnetic simulations using OOMMF code [22[. The 
equilibrium configuration of the DW in the Py layer given 
by the OOMMF simulations is shown in Fig. 1(c). 

We excite motion of the DW trapped in the IR state 
by applying an alternating current with rms amplitude 
lac between the top and the bottom leads of the device. 
This current applies ac ST to the DW and induces oscil- 
lations of the DW position along the nanowire. The DW 
oscillations give rise to spin valve resistance oscillations, 
SRac, due to the GMR effect, and a rectified voltage, 
Vdc ~ lacSRac, is generated by the device. Measure- 
ments of Vdc give information on the current-driven DW 
velocity. In our measurements, we sweep the frequency of 
the ac current, /, in the 0.5 - 10 GHz range and measure 
Vdc as a function of frequency [H , 24 1 . The amplitude of 
the DW oscillations reaches a maximum at the resonance 
frequency, /o, of the DW determined by the strength of 
the pinning potential, and a resonance peak is observed 
in the measured response curve Vddf)- 

Figure 2(d) shows Vddf) response curves measured in 
the middle of the field range where the IR state is stable 
{H = -50 Oe) with the DW trapped in the pinning po- 
tential and Bit H = 1 kOe with uniform magnetization of 
both Py and CosoFeso layers. The -50 Oe response curve 
shows a resonance at 4 GHz due to the DW oscillations in 
the pinning potential. We also observe a low-frequency 
signal, which we attribute to the DW creep in inhomo- 
geneities of the pinning potential [l^, Slj . 

The dependence of the DW velocity on current is given 
by the measurements of Vddf) as a function of lac- 
The maximum DW velocity Vmax = 27r/oAz(/o), where 
Az{fo) is the amplitude of DW oscillations at resonance, 
is obtained from the data in Fig. 2(d) (23 |: 
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where Lc=500 nm is the width of the top lead and 
Ai?=37 mfi is the resistance difference between the par- 
allel and antiparallel states of the spin valve. 

We compare our data to a theoretical description 
of DW dynamics in the rigid DW approximation fisj . 
The motion of magnetization under the influence of ST 
is described by the Landau-Lifshitz-Gilbert-Slonczewski 
equation [li, [ll, [ll] : 
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where m is the unit vector in the direction of magnetiza- 
tion. He is the effective field, a is the Gilbert damp- 
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FIG. 2: (a) Magnetoresistance hysteresis loop of the CosoFeso 
/Cu/Py nanowire spin valve for magnetic field applied paral- 
lel to the nanowire axis. The Py layer undergoes magnetiza- 
tion reversal in the -|-200 to -300 Oe field range whereas the 
CosoFeso layer switches abruptly at -900 Oe. Major (b) and 
minor (c) hysteresis loops of the Py layer demonstrate that a 
DW (intermediate resistance state) can be trapped near the 
middle of the top lead in a field range from -250 Oe to -1-100 
Oe. (d) Vddf) response curve measured at -50 Oe with the 
DW trapped under the top lead and at 1 kOe with uniformly 
magnetized Py layer. A resonance at 4 GHz is observed when 
the DW is present under the top lead. The solid line shows a 
response curve fit calculated from Eq. Q with an added con- 
stant plus 1/f background. This fit assuming zero field-like 
torque (^ = 0) gives the best fit parameters P=0.65, A=2.0 
and a=0.09. (e) Response curves calculated with the same 
parameters as in (d) except for non-zero values of the field- 
like torque (^ >0). (f) The response curves calculated from 
Eq.(|4]) and from micromagnetic simulations at Jac=0.5 mA 
are nearly identical demonstrating the validity of the rigid 
DW approximation for small-amplitude DW oscillations. 
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T^TYT is ST coefficient. 



'J ~ 2deM^ 

26| . A is the ST asymmetry pa- 



ing constant of Py, 

T - 2A^ 
^ ~ (A^-|-l) + (A^-l)(m-p) 

rameter, J is the current density, P is the spin polariza- 
tion of the current, d is the Py layer thickness, e is the 
electron charge, 7 is the gyromagnetic ratio, p is the unit 
vector in the direction of current polarization and bj = 
^fflj is the field- like torque coefficient. The dimensionless 
parameter ^ describes the ratio of the field-like torque to 
ST [27[ . The spatial profile of magnetization in the DW 
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is given by [14| : 



mz = cos {d{z)) = tanh ( ^° ) (3a) 



TO^: = sm{9{z)) cos((/)) = sech J cos(0) (3b) 

Here we use a Cartesian system of coordinates with the 
z-axis paraUel to the nanowire, the x-axis in the plane 
of the Py film perpendicular to the nanowire, and the 
y-axis along the Py film normal as shown in Fig. 1. In 
Eq-dSl), zq is the coordinate of the DW center, A is the DW 
width (A= 53 nm is given by OOMMF simulations) and 
(j) is the out-of-plane tilt angle of magnetization at the 
DW center. In the rigid DW approximation (A = const), 
Eq. © can be rewritten in terms of the DW collective 
coordinates, zq and </> [l5, 18, 28 1: 



^ _ - 7a,Ci(A) = sin(20) (4a) 

-ay - + ihj - 7ajC2(A) t&n{4>) = 7-^^o (4b) 
where K= 2ttM'^ — K± is the easy-plane shape anisotropy 

dH 

constant, is the curvature of the pinning poten- 

tial, Kj_ is the easy-axis perpendicular anisotropy con- 

stant Ci(A) = and C2(A) = dc- 

scribe the effects of ST asymmetry [28j. The value of 
K± w 0.4(27rMg ) is determined from the measurements 
of resistance versus magnetic field applied perpendicular 
to the sample plane. Assuming a parabolic pinning po- 
tential 



( dz 



const j and 0^1, Eq.Q can be solved 
analytically in the absence of spin torque and damping. 
This solution describes free oscillations of the DW cen- 
ter, zo = sin(27r/ot), in the the pinning potential. 
The resonance frequency of the DW, fo, is determined 
by the pinning potential curvature [lil . [30| : 



dz 



(5) 



The solid line in Fig. 2(d) shows a fit of a theoreti- 
cal response curve Vddf) to the experimental data. The 
theoretical response curve is found by numerically solving 
Eq.® and adding a constant plus 1// noise background. 
The fitting parameters used in the fitting procedure are 
a, P and A. The fit assuming zero field-like torque 
(^ = 0) gives ^=0.65, A=2.0 and a=0.09. This value of 
the damping parameter significantly exceeds the typical 
Py value a wO.Ol. We directly measure damping of the 
lowest-energy spin wave mode of the Py nanowire in the 
state of saturated magnetization using a ST ferromag- 
netic resonance technique [23, 24, 31| and find a=0.09. 
This shows that the high value of a is an intrinsic prop- 
erty of the ultra-thin Py film and does not originate from 
emission of spin waves by the moving DW. The observed 
large a can be explained by diffusion of Cu into the Py 
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FIG. 3: (a) Vddf) response curves measured at several mag- 
netic fields, H±, applied in the plane of the Py layer per- 
pendicular to the nanowire axis, (b) Frequency of the DW 
resonance as a function of H±. Solid line is obtained from 
micromagnetic simulations. The hysteresis in the DW reso- 
nance frequency as a function H± is due to the field-induced 
reversal of the DW chirality at -ffx«350 Oe. 



film known to increase damping [32| and by spin pump- 
ing important for ultra-thin ferromagnetic films [33| . The 
theoretical fit curve in Fig. 2(e) shows deviations from 
the experimental data at frequencies above the DW reso- 
nance. We attribute these deviations to excitation of spin 
waves in the Py nanowire. We calculate the bottom of 
the dispersion relation of the lowest-frequency spin wave 
mode in the Py nanowire to be at 5.5 GHz [3l| in agree- 
ment with the data in Fig. 2(d). 

Figure 2(e) shows theoretical response curves for a non- 
zero field-like term (^ > 0) in the equations of motion. 
Our calculations show that a non-zero field-like torque 
adds an antisymmetric Lorcntzian component to the re- 
sponse curve and decreases the quality of the theoretical 
fit to the experimental data. Fig. 2(e) illustrates that 
theoretical Vddf) curves with ^ > 0.1 show substantial 
deviations from our experimental data and thus our mea- 
surements set an upper limit on the value of field-like 
torque in our spin valve structure. 

Wc also compare our experimental data to results of 
micromagnetic simulations. The squares in Fig. 2(f) 
give the response curve obtained from OOMMF simula- 
tions. In these simulations, a spatially-nonuniform field 

OH 

Hp = — gf'Z that mimics the 4 GHz pinning potential 
is applied parallel to the Py nanowire and ST from lac 
~ 0.5 mA is applied to drive oscillations of the DW in 
the nanowire. The solid line in Fig. 2(f) obtained by 
numerically solving Eq.@ is in a good agreement with 
the micromagnetic simulation results. 

To further confirm that the observed resonance in 
Vddf) is due to DW osciUations and not due to spin 
wave excitations, we measure the resonance frequency as 
a function of magnetic field, H±, applied in the plane 
of the Py film perpendicular to the nanowire axis. This 
magnetic field modifies the curvature of the pinning po- 
tential, thereby changing the DW resonance frequency 
[loj . Fig. 3(a) shows Vddf) rneasured at several values 
of H±, and Fig. 3(b) shows the dependence of fo on 
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FIG. 4: (a) Vdc(f) response curves measured at several values 
of the excitation current, lac- (b) Maximum speed of the 
oscillating DW in the pinning potential versus the rms ac 
drive current density, Jac- Squares are the experimental data, 
dashed line is the rigid DW approximation result given by Eq. 
(|4]), solid line is given by micromagnetic simulations. 



H±. The resonance frequency shifts with H±, and the 
frequency shift is odd in H± . This odd frequency shift is 
a clear signature of the DW resonance for which the sign 
of the shift is determined by the DW chirality • A fre- 
quency shift symmetric with respect to H± is expected 
from spin wave resonances in this geometry and thus the 
data in Fig. (3) exclude the possibility of spin wave origin 
of the observed resonance. The hysteretic jumps of /o at 
H± « ±350 Oe correspond to switching of the DW chi- 
rality induced by H± . The dependence of /o on H± given 
by micromagnetic simulations is shown in Fig. 3(b). It 
is in good agreement with the experimental data. 

The dependence of the DW velocity on current is given 
by measurements of Vdc{f) as a function of lac shown in 
Fig. 4(a). Fig. 4(b) shows the maximum instantaneous 
DW velocity, Vmax, versus current density, Jac, calculated 
from the data in Fig. 4(a) using Eq.([T]). We observe 
velocities as high as 800 m/s for Jac=9xl0^ A/cm^. The 
amplitude of the DW oscillations at Vmax — 800 m/s 
is 32 nm while the smallest DW oscillation amplitude 
detectable with our technique is 1 nm. 

Because the DW mass /i ~ Mg and the ST coeffi- 
cient Qj ~ i/Ms, we expect the DW velocity to scale 
as . The DW velocity at resonance is also expected 

to be inversely proportional to damping Vmax ~ 
Therefore, for standard Py (Ni8oFe2o) material parame- 
ters (Ms=800 emu/cm"^, a = 0.01), the expected max- 
imum oscillating DW velocity in the CPP geometry is 
approximately 2000 m/s at the current density of 10^ 
A/cm^. DW velocities as high as 325 m/s at a current 
density of 5 x 10^ A/cm^ were observed in the current- 
in-plane geometry for samples with the standard Py ma- 
terial parameters (33 |. 

The rigid DW approximation predicts a linear increase 
of the DW velocity with Jac but the experimental data in 
Fig. 4(b) show clear deviations from the linear behavior 
for Jac > 1.5x10^ A/cm^. The origin of these deviations 
is revealed by micromagnetic simulations. The results of 
these simulations for Vmax{Jac) are shown in Fig. 4(b). 



The sublinear dependence of the DW velocity on current 
similar to that observed in the experiment arises from 
the breakdown of the rigid DW approximation. The mi- 
cromagnetic simulations show significant deformations of 
the DW shape induced by the pinning potential at large 
amplitudes of DW oscillations. 

In conclusion, we make measurements of DW dynam- 
ics in Py nanowires excited by spin-polarized current ap- 
plied perpendicular to the nanowire. In this geometry, 
we observe DW velocities as high as 800 m/s at current 
densities of 9x10^ A/cm^. The high DW velocities ex- 
cited by current applied perpendicular to ferromagnetic 
nanowires are promising for nonvolatile memory and mi- 
crowave signal processing applications 35|. This work 
was supported by the NSF Grants DMR-0748810 and 
ECCS-0701458, and by the Nanoclcctronics Research Ini- 
tiative through the Western Institute of Nanoelectronics. 
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